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I present the exact analytical expression for the activation energy as a function of externally applied 
magnetic fields for a single-domain magnetic particle with uniaxial anisotropy (Stoner-Wohlfahrt 
I model) , and investigate the scaling behavior of the activation energy close to the switching boundary. 
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^ ^ A lot of effort has been spent over the last few years to understand the magnetization reversal of small magnetic 

I particlesii^'^'^'^'^'i. At sufficiently low temperatures, macroscopic quantum tunneling has been observed^iiSiii, while 
for higher temperatures thermally activated behavior may switch the magnetization of the particle^^. For the anlaysis 
of the experiments one needs to know the activation energy. Moreover, with the advancing developement of integrated 
magnetoresistive memory devices (MRAM) the dependence of the energy barrier on the magnetic field has become 
^ ' of crucial technological importance as well. In a typical MRAM array, magnetic memory cells are written by a 
I . coincident field technique, where both a selected bitline (BL) and a selected wordline (WL) create magnetic fields, 
^ ' the sum of which are strong enough to switch the memory cell, whereas the fields from either BL or WL alone are 
not sufficient to switch the cell. Nevertheless, these fields do destabilize the non-selected cells to some extent, i.e. 
they reduce the energy barrier against thermally activated switching. Also, even for the selected cells the switching at 
finite temperatures happens before the actual zero temperature boundary of stability is reached, again due to thermal 
activation during the finite duration of a write puls o^^i^^ . To estimate the life time of the information in the memory, 
one needs to know the dependence of the energy barrier on the applied fields as precisely as possible, as the energy 
, barrier enters the switching rate exponentially. 
(~| ' In the study of the switching behavior of small size magnetic particles, the Stoner Wohlfarth model plays a central 
Q [ rolei^. It describes a single domain particle with uni-axial anisotropy in an external magnetic field. The single domain 
^ approximation greatly simplifies the analysis, and becomes a good approximation if the size of the particle becomes 

[ comparable to or smaller than the exchange length, which in memory elements etched out of a thin magnetic film is 
I . typically of the order of lOOnm. The activation energy in the Stoner Wohlfarth model can be calculated trivially if 
^ ■ the external field is aligned either parallel to the preferred axis or perpendicular to it. However, so far no analytical 
solution has been known for the general case of the external field pointing in an arbitrary directioni^. Given the 
crucial importance of the field dependence of the activation energy, an exact analytical solution will be provided in 
the present paper. 
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O i II. THE STONER WOHLFARTH MODEL 

o ■ 

• • . The energy of a uniformly magnetized particle with uniaxial symmetry, characterized by the anistropy energy 
. 5^ ' density K and saturation magnetization depends on its magnetization via the angle between the magnetization 
, and the preferred axis, 
;_i ' 

; E{^) = KV sin^ d - VMsH^ cos d - VMsHy sin , (1) 

where and Hy are the magnetic field components parallel and perpendicular to the preferred axis in a plane 
containing the magnetization and the preferred axis, respectively, and V denotes the volume of the samploi^. In the 
following dimensionless variables will be used by writing the energy in terms of 2KV , e{{}) = £{■&)/ {2KV), and the 
magnetic field components in terms of the switching field He — 2K/Ms as hx = H^/ hy = Hy/Hc, such that 

6(1?) = - sin^ 1} — cos d — hy sin ?9 . (2) 
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For vanishing magnetic fields, this model has a bistable ground state, whereas for very large fields {hx ^ 1 or 
hy 3> 1) the first term can be neglected, and there is only one minimum in the interval — tt < i5 < tt, such that the 
magnetization tends to align parallel to the applied field. Thus, for increasing field, one of the original minima has 
to disappear, and the fields where this happens mark the boundary of bistability. These fields are easily obtained by 
setting both the first and second derivative of ([2]) to zero and eliminating -d, whereupon the famous Stoner-Wohlfahrt 
astroid hx — cos^ hy — sin'^ is obtainedi^. 



III. ACTIVATION ENERGY 



The activation energy can in principle be calculated in a straight forward manner by finding the two minima and 
the two maxima of the energy in the bistability range, determining the metastable of the two minima and the energy 
barrier to the smaller of the two maxima^^. In the case oi hx = or hy = this is straight forward: For hy — 0, 



de 
d-d 



— cos d sin + hx sin d - 



hy COS — 



(3) 



leads to d — 0, ■& — TT, or cosz9 — —hx- The solution = is metastable for —l<hx< 0, stable for hx > 0, and 
unstable for hx < —1. Correspondingly, i) — tt is metastable for < hx < 1, stable for hx < 0, and unstable for 
hx > 1- The solution cos-d = —hx leads to two maxima which become complex for \hx\ > 1. Thus, the energy barrier 
for switching from = tt to = is given by Ea = e(arccos(— /la,)) — e(7r) = (1 — hxY /2. Similarly, for hx — 
one finds the energy barrier Ea = (1 — hy)^ jl. Therefore, if one of the two field components vanishes, the activation 

energy depends quadratically on the distance from the stability boundary. 

In the general case, where neither h X nor hy vanish, one may substitute sini9 — li, cosi9 — zt^/l into (|3|l . 
Squaring the equation leads to 



hy) — 2uhy 



0. 



(4) 



In order to calculate the energy barrier, one needs to find the roots of this 4th order polynomial, which makes the 
analysis much more cumbersome than for hxhy = 0. Still, the roots of a 4th order polynomial can be obtained 
analytically. The solutions -di are conveniently written in terms of the functions 
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Eq.(|3]) has four solutions for u, but the ambiguity in the sign of cos ■& in terms of u leads at this point to eight solutions 
for {}. They differ by three signs /i, v, a in various places, and are given by 



a arccos 



hx U v 
T+2^^+2\ 
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(9) 



The signs will be chosen according to the binary decomposition (/z, a) ~ (— , — , — ) for i^i, (— , — , +) for i92, (— , +, —) 
for ??3, . . ., (+, +, +) for i?8. None of these solutions solves (O for all hx,hy. Rather, each of them solves the equation 
only in two quadrants. However, it is possible to construct four uniform solutions out of the eight partially valid ones, 
which are continuous (after the identification of tt with —tt) and solve ([3]) for all values of hx,hy. The four uniform 
solutions are 



k = e{hy)d2 + e[-hy)di , 

h = e{hxhy)d^ + e{-hxhy)dz , 

^3 = e{hxhy)'di + e{-hxhy)dQ, 

= e{hy)'d^ + e{-hy)'&T , 



(10) 

(11) 

(12) 

(13) 
(14) 
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where 9{x) = 1 for a; > and zero elsewhere denotes the Heaviside function. Calculating the second derivative one 
finds that the first uniform solution is a minimum for hj. < 0. It disappears as real solution for > outside the 
astroid. Thus, this is the metastable state for > 0. On the other hand, 1)4 remains a minimum even outside the 
astroid for > 0, but disappears as real solution outside the astroid for < 0. The uniform solutions 1^2 and i?3 
are both maxima inside the astroid, and become complex outside one half of the astroid (for hy < and hy > 0, 
respectively ??). One easily convinces oneself that for hy > 0, is the relevant maximum (e(?93) < e{'&2)), whereas 
for hy < 0, escape over the maximum at i?2 is dominant (e(i?3) > e(i?2))- The activation energy out of the metastable 
state is therefore given by 

Ea = e{^3) ^ e0i) ioT hy > , (15) 
Ea = e(i?2) - e(z9i) for hy < , (16) 

FiglD shows a plot of the activation energy as function of h^ and hy. 




FIG. 1: Activation energy as function of hx and hy, out of the state stable for hx < 0. Outside the astroid there is only one 
stable and no metastable state, and the activation energy is therefore undefined. 



For hy = both maxima lead to the same activation energy, and Ea is an even function of hy, as it should be. In 
the following the attention will therefore focus solely on the case hy > 0, and to the case where the initial state i? = tt 
is metastable, hx > 0, i.e. the first quadrant. In principle, the particle might be excited also out of the stable state 
and end up (for a finite time) in the metastable state, but this process is of much less importance, as the activation 
energy is much larger and it enters exponentially in the thermal switching rate. 

IV. SCALING 

The scaling of the activation energy as function of the distance 1 — h from the astroid boundary is of particular 
interest. Here, h is defined by hx = hcos'^ ^, hy = hsm^ ^ such that h = 1 always corresponds to the astroid boundary. 
It has been shown beforsil that the scaling must be of the form 



Ea = cy.iOil - hf/' + C2i0il - h)' + c,/,{0{l - hf/' + 



(17) 
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All coefficients besides C2{£,) vanish for ^ = 0, and but it has been shown numericallyi^ and theoreticallj^ that 
already for small values of ^ the coefficient C3/2 dominates the scaling. This is confirmed and made more precise by 
using the exact solution. Fig. [2] shows InE'^ as function of ln(l — h). The plot reveals power law behavior to a good 
approximation even rather far away from the astroid boundary, Ea oc (1 — h)"" with an exponent a = 2 for ^ = 0, and 
an exponent close to 3/2 for larger values of ^. 




FIG. 2: Scaling of the activation energy Ea as function of the relative distance h from tlie astroid boundary [hx — /icos'^^, 
hy = ftsin"^ ^). The curves show decreasing slopes for values of a: = ^/tt = (red), 1/24, (pink) 2/24, . . . 1/4 (blue). The dashed 
line corresponds to a power law with exponent 3/2, the dot-dashed line to an exonent 2. 



According to the numerical evaluation of (fTS]) the exponent is symmetric with respect to ^ = 7r/4. The exact 
value of the exponent depends on the fitting range and on whether a quadratic term is included in the fit of In Ea as 
function of ln(l — h). Fig. [3] shows the fitted exponent a as a function of x = ^/tt for three different fitting ranges, 
from h = 0.932 to h — 0.99999853, assuming a pure power law, XuEa — aln(l — h). The observed dependence of the 
exponent on x is very similar to what was previously calculated numericallj^. In particular, the exponent appears 
to become even slightly smaller than 3/2 for values of ^ close to 7r/4. 

However, there is a substantial non-linear part in the scaling behavior, as becomes obvious when fitting to ItiEa = 
aln(l — h) + 61n^(l — h). The extracted linear part is plotted in Figd) For small values of i^, the exponent can now 
be substantially larger than 2. 



V. SUMMARY 

I have derived the exact analytical expression of the activation energy for single domain switching of small magnetic 
particles in arbitrary magnetic fields (Stoner-Wohlfarth model). The activation energy scales approximately like a 
power law as a function of the distance of the switching boundary (astroid) up to distances of order unity, but also 
contains a substantial non-power law term. 
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FIG. 3: Scaling exponent a extracted from a fit of IxiEa to aln(l — h), using the exact expression for Ea, eq. (|15p . The three 
curves correspond to fitting ranges h = 0.93187079 to h = 0.996837 (top curve — blue), h = 0.99538 to h = 0.999978 (middle 
curve — green), and h = 0.9996837 to h = 0.999998532 (bottom curve — red). For all calculations of the activation energy, 80 
digits precision was used. 




FIG. 4: Scaling exponent a extracted from a fit of In Ea to a ln(l — ft) + 6 ln^(l — h). The three curves correspond to the same 
fitting ranges used in FigO 
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